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Abstract. We find a solution to Einstein field equations for a regular toroidal lattice 
of size L with equal masses M at the centre of each cell; this solution is exact at order 
M/L. Such a solution is convenient to study the dynamics of an assembly of galaxy-like 
objects. We find that the solution is expanding (or contracting) in exactly the same way 
as the solution of a Friedman-Lemaitre-Robertson- Walker Universe with dust having 
the same average density as our model. This points towards the absence of backreaction 
in a Universe filled with an infinite number of objects, and this validates the fluid 
approximation, as far as dynamics is concerned, and at the level of approximation 
considered in this work. 



PACS numbers: 04.20.-q, 04.20.-Cv, 98.80.- Jk 
1. Introduction 

In modern cosmology, the large scale structure of the Universe is usually described by a 
Friedmann-Lemaitre-Robertson- Walker (FLRW) geometry, associated, via the Einstein 
field equations, with a perfectly homogeneous and isotropic distribution of matter. The 
FLRW geometry supplemented by perturbation theory proves to be an excellent model 
for the Universe on large scales, considering the wide range of independent observations 
that it can fit. Nevertheless, this ability to fit most, if not all, the observations available 
today comes at a cost: the unavoidable introduction of dark matter and dark energy: 
two elusive energy components that actually dominate the dynamics of the late time 
Universe on all scales. The introduction, by hand, of these new, exotic, components 
thus raises the concern that something may be wrong in the more fundamental way we 
describe the late-time Universe, i.e. it inevitably leads one to wonder if our theory of 
gravitation and/or our assumptions on the geometry of the Universe might be oversim- 
plified. The first possibility advocates for a modification of the laws of gravity on large 
scale (see, for example [1] for a recent review). Weaker assumptions on the geometry of 
the Universe have also been explored, mainly in the context of Lemaitre-Tolman-Bondi 
models; see, for example |2l El HI O [6l [7] , but also [H [9] for criticisms of this possibility. 
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Of course, the actual Universe, at least in its later stage of evolution, is not perfectly 
homogeneous and isotropic, because of the gravitational instability that causes struc- 
ture to form out of tiny initial fluctuations in the matter distribution, forming bound 
objects such as galaxies and clusters of galaxies. Therefore, the FLRW geometry is usu- 
ally constructed by representing, in the field equations, the complicated, inhomogeneous 
energy-momentum content by a smoothed energy-momentum tensor involving averaged, 
non local, energy densities and pressures. This is the so-called fluid approximation: on 
large scales, the energy content of the Universe is well modelled by one or several ho- 
mogeneous and isotropic perfect fluids. Despite the fact that it is widely used, it is 
not very often stated explicitly, let alone constructed (see for a recent discussion). 
Moreover, it is now well-known that averaging cosmological quantities leads to the in- 
troduction of backreaction terms [HI |12l UHl HI] that could modify the kinematics of the 
model on large scales, even though the amplitude of such corrections is not presently de- 
termined and is still a matter of debate (see [15] for a critical discussion on these issues). 

The purpose of this paper is not directly to discuss the backreaction issue. Rather, 
we would like to propose a solution of Einstein field equations corresponding to a set of 
point masses representing objects such as galaxies, in gravitational interaction with each 
other in an otherwise empty spacetime. This will allow us to discuss the kinematics of 
the collection of point masses and compare it to the kinematics of an analogue FLRW 
model with a fluid corresponding to the one emerging from the distribution of masses 
when it is smoothed. This route has recently been explored by Clifton and Ferreira 
in a series of papers [161 IH] based on the Lindquist- Wheeler tessellation of spacetime 
[18] . The Lindquist- Wheeler approximation consists in pasting together a large number 
of Schwarzschild cells representing the metric around point masses. Then, the Israel 
junction conditions imply that the boundary on which the cells are glued expands fol- 
lowing the FLRW kinematics. Unfortunately, it is impossible to glue these cells together 
exactly: the Israel junction conditions are only satisfied on a subset of the boundaries, 
where they intersect, but the boundaries cannot be exactly matched, leading to an over- 
lap region and a 'no-man's land' [16]. Also, spherical symmetry in each cell is only 
approximate. Therefore, despite its elegance, this method relies on approximations that 
are difficult to control, and one would hope for symmetric situations in which an exact 
solution could be found. This was done in [19], where two antipodal masses were put on 
a 3-sphere. It was shown that an exact solution to this problem could only be found by 
introducing a cosmological constant. Recently, and independently of the current work, 
[20] have studied a more realistic solution in the context of cosmology, namely, a finite 
number of masses arranged in a regular lattice on a hypersurface having the topology 
of a 3-sphere. Nevertheless, they only studied the kinematical backreaction arising at 
a given time (the time of maximum expansion) by solving the constraints equations, 
leaving the full dynamics for further studies. A similar study, but for the fiat topology 
was also conducted numerically in [2lj . 
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In this paper, we will show that, by considering an infinite cubic lattice of size L, with 
equal masses M at the centre of each cell, one can find a solution to the Einstein field 
equations that is exact at second order in an expansion in powers of y^MjL. Such an 
approximation scheme improves on the Lindquist- Wheeler approximation, since there is 
no need for junction conditions (and hence, no 'no-man's land'), and because it retains 
a certain degree of anisotropics around each mass. The main result of the paper is thus 
the expression for the metric, Eq. f l22|) . A similar type of solution was studied in [22] . 
with periodicity along one axis. [23] also investigated such lattices, by concentrating 
on small deformations of the geometry of each cell in powers of v/c, the velocity of the 
objects described by the lattice. Our approach is very similar to this work, but instead 
of solving the perturbed junction conditions between cells, we solve directly for the full 
metric of spacetime perturbatively. 

The choice of ^J MjL as a small parameter will be motivated by the fact that a regular 
Taylor expansion of the metric in powers of Mj L cannot satisfy the field equations in a 
toroidal symmetry. Rather, one has to rely on an expansion in terms of M/L, because 
in this case, the source terms in the equations, that is of order M/L, can be reabsorbed 
by a correct re-definition of the integration constant at order ^J M/L, leading to a well- 
behaved solution. It will also be shown that, by the same process, the scalar solution 
at order M/L is fully determined by requiring that the equations be satisfied at order 
(M/L)^'^^. In Sect. [2], we present the solution obtained when the distribution of masses 
is modelled by a Dirac comb. We discuss the subtleties associated with the series expan- 
sion and present the method used to obtain the solution in details. We will see that the 
Fourier expansion of the solution leads to diverging terms in the metric in the form of 
non convergent series exhibiting ultraviolet divergencies. Therefore, there is no solution 
to the problem for a singular distribution of masses, in agreement with the results of 
[22] . In Sect. [3l we address this non- convergence problem by smearing the distribution 
of mass. This amounts to replacing the Dirac comb by a sum of approximations of 
the Dirac 5 function via peaked Gaussian functions. This will introduce a natural UV 
cut-off in the diverging parts of the metric, and this cut-off will be set properly to obtain 
a solution that describes accurately the spacetime exterior to the masses. In Sect. HI we 
discuss the kinematical characteristics of the solution via its the effective scale factor, 
and we compare this one to the analogue FLRW model with dust, where the energy 
density of the dust fluid is identical to the one derived from the distribution of point-like 
mass in our lattice. We flnd that the fluid approximation is valid for the special solution 
presented in this paper, since the scale factors derived from the exact solution and from 
the analogue FLRW model match at order (M/L)^/^. This is in agreement with the 
results of [20], since they found that the kinematical backreaction in their flnite model 
decreases as the number of masses increases: here we have an inflnite number of masses, 
so that at the limit, we would indeed expect the backreaction effect to exactly cancel. 
Nevertheless, this comparison should be made with caution, since we only have a per- 
turbative solution at order (M/L)^/^; non-linear effects at higher orders could generate 
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backreaction through mode couphngs. Finally, Sect. serves as a concluding section 
and introduces some future works. 

We use the sign convention (— , +, +, +) for the spacetime metric. Latin indices are used 
for space components, and boldface symbols denote three vectors as well as triplets of 
integers. We also use natural units G = c = 1 throughout the paper, except stated 
otherwise. 

2. A lattice Universe made of point masses 
2.1. Framework 

Physically, we are trying to describe the collective behaviour of typical astrophysical 
objects such as galaxies. This is clearly an impossible task if we do not introduce fur- 
ther assumptions on the symmetries of the problem. The idea is to obtain a model that 
is genuinely inhomogeneous, but exhibits an averaged homogeneity when considered on 
large scales. A regular lattice of identical masses certainly possesses these character- 
istics, and we will consider a cubic lattice of size L, with identical masses M at the 
centre of each cell of the lattice. Therefore, space clearly has a toroidal symmetry. The 
periodicity of the distribution of sources will be extremely valuable when trying to solve 
the field equations, since it will impose a periodic solution, and allow the use of Fourier 
series. 

Nevertheless, even formulated as such, it is extremely difficult to solve the full problem, 
because of the non-linearity of the field equations. But we can do one further assump- 
tion based on the physics we are trying to understand. If the masses on the lattice are 
to represent typical galaxies, we can choose, as our typical parameters 



where M© ~ 10^° kg is the Solar mass, and L is of order of the intergalactic distances. 
With these numbers, the natural parameter of the lattice is Rs/L ~ 10^^ ^ 1, where 
Rs = 2M is the Schwarzschild radius of the masses. Therefore, we can look for a solution 
expanded into powers of M/L; this will lead to linearised field equations that can be 
solved exactly, and we expect that the first few terms in the solution will represent the 
spacetime exterior to the masses with high accuracy. 

2.2. Definition of the perturbative expansion 

In the following we will consider a source term given by a three dimensional Dirac comb 
in comoving coordinates Too oc M'^^^^a 5^^\x — Ln), where the masses are located at 
Xn = Ln. Such a periodic distribution of masses can be written as a Fourier series with 



M ~ lO^^Mp, and L ~ 1 



Mpc, 




(1) 
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showing that the zero mode of the source term is a constant nonzero (comoving) density 
p = M/L^. We could look for a periodic solution for the metric of the standard form: 
dfiu = Vfj.f + f^f^fj where h^u would be of order M/L. However in this case, one finds, 
Fourier expanding h^^j,, that G^q^*^ mode _ 0{{M / LY) showing that no periodic metric 
of the form considered above can solve Einstein equations for such a matter source. 

This is reminiscent of an analogue situation for the Coulomb potential in infinite lattice 
of charges, i.e. crystals. In crystals the zero mode of the periodic Coulomb potential 
cannot be solved in general, and one simply drops this term by noticing that the average 
density of charges vanishes, provided the crystal is electrically neutral [211 [25] , see also 
|26j . In our case however, the global scalar charge does not vanish and the zero mode of 
Einstein equations, in a periodic solution, must be absorbed through an homogeneous 
expansion of space. Comparing to a flat and homogeneous, dust-dominated universe, we 
see that Hubble constant relates to the density of matter via H'^{t) = Svrp/S, showing 
that the corresponding FLRW metric contains a ^/p ~ a/ M/L^ term. 

As a consequence, we need instead to consider a perturbative expansion for the metric 
of the following form 

9,^ = V,u + ^ -j^h^ + -j^hf, + + O (^-^ J , (2) 

where the first order term will absorb the zero mode of the source term through an 
homogeneous expansion of space (i.e. we will select the solution for which h^i^} depends 
solely on time), whereas the anisotropic gravitational field created by the masses will 
be described in the h)j,v. part. This makes clear from the beginning that a solution 
up to the order 3/2 in M is attainable, mixing the global expansion with the local 
anisotropic gravitational field. However the order solution is out of reach since 
Einstein equations at this order involve non linearities, and in particular product of 
Fourier sums and non-trivial mixing of the Fourier modes. With the ansatz Eq. ([2]), we 
are now in position to derive both the expressions for the matter source and linearised 
field equations. 



2.3. Matter source 

The stress-energy tensor corresponding to our lattice universe reads 

Comoving coordinates require = and therefore = 5^/ a/— (?oo(xn)- Energy 
conservation imposes Vfj_T'^^ = 0, and is automatically solved to any order by choosing 
the synchronous comoving gauge goo = —1 and goi = in which we will work for now 
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on. Then one shows that the only non-vanishing component of T^^^ to order M^/^ reads 

Too = M(l-'-Mj2hi^j:5^H.-Ln), (4) 

\ i J nGZ3 

where we used Eq. ([3]), the gauge conditions, the fact that T^'-' is 0{M) to the leading 

order, and (7^1, = 7]^,, + sj^h^^l + 0{M). Fourier expanding the Dirac Comb, we finally 
get 



L3 I 2 V L ^ " / ^ V L2 

In the following we shall also work with the wave vector k = 21111/ L and write 

Metric, field equations and solution 

In the synchronous gauge comoving with the masses, we split the line element w.r.t. 
scalar, vector, and tensorial parts 

ds'^ = -dt^ + [(1 + 2C)5ij + 2dlE + 2d(^iEj) + Eij] dx'dx\ (6) 

where, according to Eq. ([2]), we write C = ^/mJLC^^^ + {M/L)C^^'> + {M/Lf/^C^^'^ 
and similarly for the other quantities. Then divergence free and tracelessness split 
accordingly into diE^^P'^ = 0, diE^^'^^^ = and E^^^' = for p = (1,2,3). Einstein 
equations are decomposed into their scalar, vector and tensor parts. 

2.4- 1- Order a/ M/ L field equations We get, for all {i,j), and where a prime denotes 
a derivative with respect to t 



Scalar part 



= AC(^) (7) 
= diC^^^' (8) 
= (E(i)" - - A - - 2C(^)"5,„ (9) 



Vector part 



= AeI^^' (10) 



= dM^" + d,M'^", (11) 



Tensor part 



Eff" - = 0. (12) 

We look for periodic functions and write, e.g. C^^^(t,x.) = C^^\t) + ^ii^oC^\t)e^^^. 
The system of equations is then solved mode-by-mode and one easily gets the most 
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general solution 

L 

^(1) = L'E^^\t) + LK{^)t + LV(x) 

where K and J are any dimensionless zero mode free periodic functions. E^^\t) and 
E^^\t) are not determined but do not appear in the expression for the metric, and 
can thus be taken to zero without loss of generality. Physically, we see on the above 
equations that the vector modes are not sourced and do not propagate, explaining why 
e[^^ = 0. On the other hand the tensor modes propagate but are not sourced, so that 
we recover the standard gravitational waves. 

Derivatives of the functions K and J appear in the metric. These terms have to be 
understood as being the most general source terms compatible with the symmetries and 
which do not contribute to Einstein field equations to this order. They are "spurious" 
source terms with no physical counterparts. Indeed it is easily shown that they act as 
source terms for order M equations, and enter these equations in a non linear way. It 
is thus physically meaningful to set them to zero (and anyway necessary to be able to 
proceed further in the perturbative expansion). The same remark applies to the grav- 
itational waves found above. First they have no clear physical origin, as they are of 
order Vm, and second, if not dismissed, they would source the tensorial field equations 
to order M in a non trivial way. 

We thus restrict ourselves to the simplest solution, namely C*-^-* = A/Lt + B, E^^^ = 0, 
E^^^ = and E^j^ = and proceed to the next order. As discussed previously, the 
solution we picked up only depends on t and essentially absorbs the zero mode of the 
source term at order M, via the constant A. This is what we show now. 



2.4-2. Order M/L field equations Let us start with the vector and tensorial equations. 
Using the previous a/M solution, they are found to be the same as Eqs. ffTOl [TT| [T2|) up 

— ^2) — (2) 

to the fact that it now applies to second order quantities E^ and E^j . Thus, the vector 
part vanishes, while we choose to dismiss again the gravitational waveiji. We therefore 
focus on the scalar equations, which are 

3^-2AC<^) = ?^5:e- (13) 

k 

I One could keep these gravitational waves, as they will only source the order field equations, which 
go beyond the present solution. However we are more interested in the gravitational field created by the 
masses and the subsequent expansion of the lattice rather than by its possible content in gravitational 
waves. 
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diC^^^' = (14) 



(15) 



We see how the order ^jMjL solution for C {C^^^ = At/L + B) now enter the order 
M/L equations, and can absorb the zero mode of the source term (to this order). Zero 
and non zero modes of Eq. (fT3|) indeed give, respectively 

A^e^ (16) 
CP'((,x) = CP'(«) + /(x), (17) 
where e = ±1 and where the function /, central to our solution, is given by 



TT ^-^„ In 



• 27rn.x , 

^e^^, (1^ 



and will be discussed below. The star in the notation means that the null triplet 
(0, 0, 0) is dismissed. Equivalently we have 

f(^\-l cos nx) cos {^py) cos {^qz) 

(n,p,g)GN3 ^ ^ 

Eq. f|T4|) is then trivially satisfied, while Eq. (|T5|) gives both 

E^'\t, x) = L2E(2)(t) + |/(x) + Lf/(x)t + LV(x), (21) 

where again E^'^\t) is left unspecified but drops from the metric, whereas now U and V 
are undetermined arbitrary (periodic, zero mode free) functions. 

2.4-3. Order (M/L)^/^ field equations We will not detail here the equations of motion 
as they are more complicated. One can easily derive them using any tensorial package, 
and solve them in a similar fashion as was done for the previous orders. This is 
straightforward although a bit tedious. Focusing on the scalar part, we have been 
able to derive the most general (scalar) solution which is displayed in the Appendix, 
see Eq. flA.ip . It is the "most general solution" provided that one sets K and J found 



in the a/M solution to zero, and dismissing the gravitational waves. In the course of 
solving order 3/2 equations, it turns out that, in analogy with the transition between 
the solution at order a/ M/L and the solution at order M/L, during which the constant 
A acquired a value in order for the equations at order M/ L to be satisfied, the constant 
a in Eq. ( l20l) has to take a specific value 

a = eBJ-. 



The integration constants B, /3, T, and w found in Eq. ( lA.ip . if not dismissed, would 



enter the effective scale factor (see Sect. H]) at post-Newtonian order. These terms have 
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thus no Friedmannian counterparts and it is natural to set them to zero. The arbitrary 
functions U, V, R and S (see Eq. flA.ll) ) affect the local dynamics of spacetime, but again, 
are essentially source terms with no physical counterparts and would be constrained by 
the higher order field equations. We thus let them to zero and obtain the simplest 
(M/L)^/^ solution for the lattice universe (putting back G and c): goo = —1, goi = 0, 
and 



9ij 



l + 2e 



[gmV'^ 



Lc2 




(22) 



3 3L 

where the function / was given in Eq. ( fT9|) above, and we remind that e = ±1. We 
will show in Sect. Hlthat e relates to the sign of the effective Hubble constant: e = 1 
corresponds to an initially expanding universe that decelerates under the effect of grav- 
ity, while e = — 1 is an initially contracting Universe that accelerates its contraction. 
This metric solves Einstein equations for the matter source given in Eq. (jl]) up to order 
{M/Lfl\ 



Clearly the solution is valid only in the limit where g^i, — rj^y <^ 1. Under the identi- 
fication of the effective Hubble constant H = ^ySn/Sy^ M/ (see Sect. Hj), the time- 
dependent terms in the above metric show that our solution is valid locally in time 
At < H-^ oc L^/L/M, while space-dependent terms are more complicated to anal- 
yse because of the non trivial form of the function /. Still, based on dimensional 
analysis at least, we also expect our solution to be valid only over a similar range 
Ax ^ H^^ oc Ly^L/M. Note that, for the numbers quoted previously, corresponding 
to a lattice of galaxy-like objects, the time span valid for the solution is of order 1 Gyr, 
which is already a good cosmological time, allowing the use of this solution for cosmo- 
logical predictions. 

The function / and its derivatives are divergent, and so is the metric. However next 
section will show how to regularize it by introducing a physically motivated cut-off. It 
is enough, here, to note that the divergence in / comes from large values of the Fourier 
modes and is therefore a UV divergence linked to the point-like nature of the sources. 
Therefore, if the sources are extended but still very small compared to L, we expect to 
obtain a similar solution. This will indeed be proven in the next section. 



3. A Lattice universe with peaked Gaussians 



The function /(x) and its derivatives are ill-defined, because the series expansion that 
defines /(x) does not converge. However, this problem comes from a UV divergence 
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of the expansion, and can be easily dealt with. One should remember that one can 
approximate the Dirac 6 function by 

6(x-nL) -e^^, (23) 

for 7] arbitrarily small, and where rj has a dimension of a length. This rj will be related to 
a UV cut-off applied to Fourier modes. Then the (total) regularized source term reads 
S{x,y,z) = S{x)S{y)S{z) with 

(24) 

riJn 

Such a source term turns out to be the convolution of a Dirac Comb with a normal 
distribution of variance ■i]'^. By virtue of the convolution theorem its Fourier transform 
is thus given by the product of the Fourier transforms of the Dirac Comb and the normal 
distribution. One then finds that the total source term can equivalently be written as 

Six)Siy)S{z) = j-,Y. e"^—"^. (25) 



73 



We also note that the source term is properly normalized. Indeed, taking advantage of 
Eq. ( 125|) . one can write 



oo 



^(x) = l + 2^e-^^cos ( j , (26) 

n=l ^ ' 

where one recognizes a Jacobi theta function. The above expression makes clear that 

1 /"^/^ 



L 



S{x)dx = 1, 

L/2 



since all the terms in the sum cancel. This applies as well to S{y) and S{z), so that the 
integral of the total source term on the fundamental cell [—L/2, L/2]^ is also unity. 

The source term given by the Fourier series in Eq. ( l25l) represents therefore a proper 
regular izat ion of the point-like lattice Universe we considered thus far. Moreover, the 
derivation of the solution is not affected by this change, and one can apply directly the 
same method as previously: the solution will remain the same, except that /(x) needs 
now to be replaced by /r,(x) given by 

2, |2 2 
7r |n| J7 

/.(X) = - E (27) 
TT ^ — ^ n r 

We see that r] acts as a UV cut-off, suppressing all the terms in the series with 
|np ^ • ^ natural cut-off is provided by the Schwarzschild radius of the masses: 

rj = 2M, in which case we would need the first 10^ terms of each sum. Of course, since 
we are interested in the exterior solution, we may not need to go that far. A galaxy is 
approximately 10^ times bigger than its Schwarzschild radius, therefore, for an exterior 
solution, we may take r] ~ lO^M, and then, we only need the first 200 terms of each 
sum. 
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4. The fluid approximation 

In order to characterize the dynamics of the model we calculate the rate of expansion 
between two masses. Let us consider two masses on the x-axis (all the other axes are 
equivalent, by symmetry), separated by a coordinate distance NL, for N an integer. 
The physical distance between the two masses is given by 

Kt)= / ^dx. (28) 







Then, expanding the square root to order (M/L)^/^ we find 

l{t) 



a{t) 



NL 



3/2 



8n GM 27rGMt'^ ine Sn f GM\ , 

because the part of the metric that has a spatial dependence is periodic, of period L, 
and where we introduced the effective scale factor a(t) of the model. The Hubble flow 
defined by H{t) = a'{t)/a(t) is then found to be (at order M^/^) 

[M IGM AnGMt , /- o,, [GMV''^ , 
H{t) = e\ —\ — -^ + 4y67r^/2£ — - ^ (30) 



3 V ^3 L3 V / 

Then we identify the initial expansion rate at t = 



/87r GM , , 

where we now choose the solution e = 1, ie. i^o > 0. This is Friedman equation for a 
flat Universe filled with dust with density p = M/L^: Hq = 87rGM/3L^. Replacing in 
the expression for H{t), one finds 

Hit) = Ho- ^H'^t + ^-H',t' + 0{H^) (32) 

and this corresponds to a flat FLRW model with equation of state w = 0. Thus, the 
model with discrete masses on a cubic lattice is identical to a FLRW model with dust, 
with the corresponding energy density. This means that one cannot distinguish the 
distribution of mass (localized or smeared in an homogeneous manner) from purely 
kinematical considerations. This is in agreement with results obtained in the Lindquist- 
Wheeler approximation [181 [E] • As mentioned in the introduction, this also agrees with 
the results of [20j where it is found that the bigger the number of masses, the smaller the 
amount of backreaction produced at a given time: here we have an infinite number of 
masses, and a zero backreaction. Of course, this result is only valid at order (M/L)^/^, 
that is at the lowest non-linear orders; mode couplings at higher orders might alter the 
picture and create some backreaction. The study of these higher orders is left for future 
works. A calculation of the averaged scale factor in the standard way it is done in 
the backreaction context [12], i.e. by using volumes averages, lead to exactly the same 
result, provided the domain respects the symmetry of the lattice (which is necessary in 
the averaging context anyway, since the domain is a Lagrangian one). 
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5. Conclusions 

In this paper we derived a new solution for a lattice cubic Universe that is exact up to 
order M/L in the lattice parameteiQ. The choice for the toroidal symmetry was some- 
how arbitrary to begin with. It is clear that the technique could be straightforwardly 
extended to cover more general "crystals" - or Bravais lattices - without any expected 
new difficultiejjj. This shall be explored in future works. 

In any event, the present solution Eq. (1221) for the cubic lattice is a significant progress 
towards the understanding of the effects of inhomogeneities on the late Universe. In- 
deed the solution covers the whole space-time directly (or at least a significant part of 
it, given by Ax <^ H~^, At <C H~^) without the need for gluing together elementary 
cells, thus eluding the intrinsic difficulties of the Lindquist- Wheeler scheme. 

There are two issues behind the late time inhomogeneities in the Universe: first, there is 
the question of their effect on the dynamics of the Universe itself. Second, it is necessary, 
especially in the era of precision cosmology, to quantify the effects of inhomogeneities 
upon the observables and in particular on the propagation of light rays; see a very in- 
teresting, recent paper on the subject p7]. In this paper we focused only on the first 
question and showed that the dynamics of the lattice Universe does justify the validity 
of the fluid approximation, in agreement with previous studies. In particular, our so- 
lution presents no backreaction at all, although it should be reminded that it is only a 
perturbative solution, valid for a limited amount of cosmic time, at a given order in the 
ratio M/L. Another limitation of this model is also clearly its regularity: in the real 
Universe, objects are structured in a more complex network made of galaxies, clusters of 
galaxies, filaments and large voids. Improving on this issue could be done in our model 
by introducing imperfections in the lattice and by studying their behaviours, as it is 
commonly done in solid state physics, through the study of defaults in crystal. This 
could be the subject of a future study. 

The question of observables in such models remains to be addressed. Is the propa- 
gation of light in the model with localized masses different from its FLRW counterpart? 
From the form of the metric, we should expect corrections at order M/L at least, due 
to the presence of the anisotropic function /. Also, and more within an astrophysical 
context, what are the effects of the inhomogeneities on the local dynamics (e.g. on 
the rotation of spiral galaxies)? Preliminary studies indeed indicate that the function 
/ found in the solution is genuinely anisotropic, even at very short scales, i.e. around 

§ Actually the solution extends to order 3/2, but we expect some of the integration constants and/or 

free functions found at that level to be fixed by the next order of perturbation. Still, Einstein equations 

are solved up to order 3/2 by Eq. (|22p . and, for instance, Eq. (|32]) shows that the effective Hubble 

constant in the lattice universe does match with its FLRW counterpart up to order 3/2. 

II At least for flat 3-geometry, to which we restricted ourselves through our choice for the Fourier 

expansions. 
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the masses themselves. Does this affect the local dynamics, and to what order of mag- 
nitude? We refer the reader to a forthcoming publication [28] on these issues for the 
lattice Universe. 
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Appendix: Most general (scalar) solution up to order Af^/^. 

The following metric at order [M/Vf'^'^ solves Einstein equations with source term given 
by Eq. (j4]). We dismissed the tensorial perturbations at any order of interest, and the 
vectorial part is vanishing. We are thus left with scalar perturbations alone, and the 
metric is 



fl'oo 






= 


9ij 


= ^ij 



2GM 



+2 



Lc2 




1 + 2 



/3 + /(x)+.5,/y^ + ^ 



Lc2 J 

2tiB cH^ 
3 L2 



5/(x) + T + ^f/(x) + ^ f + 2eJ^fi^) 



L 



27T6 8ncH 



9 




.3+3 



3 L3 



2GM^2 fcH 
-of, 



2+2 



Lc^ \ 2 

3/2 



■/(x) + LctU{x) + lV(x; 



+e 



27r cH^ 



d: 



L^S{^) + Lcti?(x) + 



cH^ 




2e\ —U{-K)-Bf\-K) 




/(x) 



(A.l) 



The quantities 5, /3, T, and w in Eq. (lA.ip are constants of integration, undetermined 
at order [M/Vf^'^. One would need to go to higher orders to determine them. The 
functions f/, V, R, S are arbitrary periodic functions with a vanishing zero-mode. 
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